Astronomy & Astrophysics manuscript no. 15072ap 
July 12, 2010 



©ESO2010 



Herschel/HIFI measurements of the ortho/para ratio in water 
towards Sagittarius B2(M) and W31C^ 

D. C. Lisi, T. G. Phillips!, p p Goldsmithi^ D. A. Neufeld\ E. Herbst'^^ q Comito^ P. Schilke^i^ H. S. P Muller'^^ 

E. A. Bergin^, M. Gerin'^, T. A. Bell', M. Emprechtinger\ J. H. Black^", G. A. Blake^ F. Boulanger^^ S. Cabrit^^^ 

E. Caux^'^ C. Ceccarelli^ J. Cernicharo', A. Coutens'^^ N. R. Crockett^, F. Danief ''^, E. Dartois^^ M. De Luca'', 

M.-L. Dubernet^''''^ P. Encrenaz^, E. Falgarone^, T. R. Geballe^^, B. Godard^, T. F. Giesen'^, J. R. Goicoechea^, 

C. Gryi^ H. Gupta", P Hennebelle^ P Hily-Blant^, R. Kolos^^ J. Krelowski^^, C. Joblin^^, D. Johnstone'^ 

M. Kazmierczak^'', S. D. Lord'^, S. Maret^, P. G. Martin^'', J. Martin-Pintado'^, G. J. Melnickl^ K. M. Menten^ 

R. Monje\ B. Mookerjea^^, P. Morris'^, J. A. Murphy'^, V. Ossenkopf'^-^^, L. Pagani^^, J. C. Pearson'^, M. Perault^, 

C. Persson^^, R. Plume^', S.-L. Qin'^, M. Salez^^^ 3 Schlemmer'^^ M Schmidt^*', P Sonnentrucker^ J. Stutzki'^, 

D. Teyssier^i, N. Trappe'^ F F S. van der Tak^", C. Vastel^-S, S. Wang^, H. W. Yorke'^ S. Yu", J. Zmuidzinas^ 

A. Boogert^'', N. Erickson^^, A. Karpov^ J. Kooi\ F W. Maiwald'^ R. Schieder'^, P ZaaP 

(Affiliations can he found after the references) 
Preprint online version: July 12, 2010 

ABSTRACT 

We present Herschel/HIFI observations of the fundamental rotational transitions of ortho- and para-H^^O and H^^^O in absorption towards 
Sagittarius B2(M) and W31C. The ortho/para ratio in water in the foreground clouds on the line of sight towards these bright continuum sources is 
generally consistent with the statistical high-temperature ratio of 3, within the observational uncertainties. However, somewhat unexpectedly, we 
derive a low ortho/para ratio of 2.35 ± 0.35, corresponding to a spin temperature of ~27 K, towards Sagittarius B2(M) at velocities of the expanding 
molecular ring. Water molecules in this region appear to have formed with, or relaxed to, an ortho/para ratio close to the value corresponding to 
the local temperature of the gas and dust. 
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1. Introduction 

Water molecules play an essential role in the physics and chem- 
istry of the dense interstellar medium (ISM). Water is one of 
the main reservoirs of oxygen, and as an important coolant of 
dense gas it strongly affects its star formation properties. For a 
molecule with two hydrogen atoms, such as water, the ortho/para 
ratio is temperature dependent and, in principle, the temperature 
' of the medium in which the proton spin state populations last 
, equilibrated can be deduced from the observations. Water is an 
asymmetric top molecule, with energy levels labeled by the set 
of quantum numbers J, K^i, K+i, where K^i, K+i are the limit- 
ing prolate and oblate symmetric top quantum numbers. Levels 
' with K-i + K+i even are para, and those with K^i + K+i odd 
' are ortho. There are no fast radiative transitions between the two 
species of water, ortho and para. However, a collision resulting 
in a proton exchange, with a proton or H3 ion, can cause an 
ortho-para conversion in a water molecule. Such conversion can 
also result from interaction between a water molecule and the 
grain surface on which it is adsorbed. The lowest ortho state of 
water is 34.2 K above the para ground state. The ortho/para ra- 
tio in water is 3 in the high temperature limit. Departures from 
this limiting value are seen for nuclear spin temperatures below 
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~50 K, with the rat io dropping to 1.5 at ~18 K (see Fig. 4 of 
lMummaetaLlll987h . 

The ground-state rotational transition of o-water at 557 GHz 
was studied extensively i n the ISM by the SWAS a nd Odin satel- 
htes dMehiick & Ber gin',l20 05|:lHialmarsonll2004l) . However, the 
HIFI instrum ent ( de Graau\\|^ 20101) aboard the Herschel Space 
Observatory (Pilbratt et al., 2010) allows for the first time het- 
erodyne studies of the fundamental p-water line at 1113 GHz. 
Measurements of the ortho/para ratio in water, by means of ab- 
sorption spectroscopy in cold foreground clouds on sightlines to- 
wards bright submillimeter continuum sources, can provide key 
insights into the thermal history of the gas. The ortho/para ratio 
in water has been measured in several solar system comets and 
the derived nuclear spin temperatures are typica lly of order 30 K 
(see lCrovisier et alJll997HKawakita et~ani2004 . 

2. Observations 

HIFI observations of the ortho and para Hj^O and Hj^O towards 
Sagittarius B2(M) and W31C presented here were carried out 
between 2010 March 1 and March 5, using the dual beam switch 
(DBS) observing mode, as part of guaranteed time key pro- 
grams HEXOS: Herschel/HIFI Observations of Extraordinary 
sources: The Orion and Sagittarius B2 Starforming Regions and 
PRISMAS: Probing Interstellar Molecules with Absorption Line 
Studies. The source coordinates are: ffy2000 = 17''47™20.35^ 
and (5y2000 = -28°23'03.0" for Sagittarius B2(M) and ayiooo = 
18M0'"28.7-^ and 5y2000 = -19°55'50.0" for W31C. The DBS 
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reference beams lie approximately 3' east and west (i.e. perpen- 
dicular to the roughly north-south elongation of Sagittarius B2). 
We used the HIFI Wide Band Spectrometer (WBS) providing a 
spectral resolution of 1.1 MHz (~0.6 kms ' at 557 GHz) over 
a 4 GHz IF bandwidth. The spectra presented here are equally- 
weighted averages of the H and V polarizations, reduced using 
HIPE (Ott 2010) with pipeline version 2.6. The resulting Level 2 
double sideband (DSB) spectra were exported to the FITS for- 
mat for a subsequent data reduction and analysis using the IRAM 
GILDAS package. 

The band la, lb and 4b spectral scans of Sagittarius B2(M) 
consist of DSB spectra with a redundancy of 8, which gives ob- 
servations of a specific lower or upper sideband frequency with 
8 different settings of the local oscillator (LO). This observing 
mode allows for the deconvolution and isolation of a single side- 
band (SSB) specti-um (Comito & Schilke, 2002). We appHed the 
standard deconvolution routine within CLASS. The observations 
of water Unes in W31C were obtained using the DBS single 
point observing mode with 3 shifted LO settings that were av- 
eraged to produce the final spectra. The HIFI beam size at 557 
and 1113 GHz is 38" and 21", respectively, with a main beam 
efiiciency equal to ~0.68. 

HIFI spectra towards strong continuum sources are af- 
fected by standing waves. For sources with narrow lines, these 
standing waves can be removed efficiently in the pipeline us- 
ing the FitHifiFringe task in HIPE. However, in the case of 
Sagittarius B2(M), the period of standing waves is comparable 
to the line widths in the band 1 data. Given the high density of 
lines in this source in the low-frequency HIFI bands, the stand- 
ing wave removal can lead to artifacts in the resulting spectra. 
We, therefore, have not attempted to "defringe" these data. 

3. Results 

Figure 1 shows spectra of the ground-state rotational transi- 
tions of o- and p-Hj''© and H2'^0 towards Sagittarius B2(M). 
Lines of Hj''© are completely saturated over a wide range of 
velocities. Absorption in the envelope of Sagittarius B2 is seen 
around 62 kms"', while absorption at velocities between -140 
and 40 km s"' is due to the foreground gas. 

The o-Hj^O transition can be observed in HIFI bands la and 
lb. While the two observations give the same line shape, the 
depth of the absorption is a few percent above the zero contin- 
uum level in band 1 a and a few percent below the zero continuum 
level in band lb (Fig. 1, upper panel). This effect is due to small 
imbalances in the mixer sideband ratios. The sideband decon- 
volution can, in principle, fit the sideband gains. However, the 
H2O line is near the band edge and can only observed in the up- 
per sideband in band la and only in the lower sideband in band 
lb. The resulting determination of the sideband gains is thus not 
completely accurate. In the following analysis, we corrected for 
this instrumental effect by assuming that the deepest absorption 
in both Hj^O lines defines the zero continuum level. 

Sagittarius B2(M) has a rich emission spectrum, which adds 
to the dust continuum, and against which the foreground clouds 
can absorb. Some of the strongest and most numerous emission 
lines seen in the spectrum are those of SO2, which we have mod- 
elled assuming LTEQThe resulting model, shown as green lines 
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Fig. 1. (Upper) Spectra of the ground-state transitions of o-Hj^O 
(band lb) and H^^O (band la) towards Sagittarius B2(M) (black 
and blue histograms, respectively). (Lower) CoiTesponding 
spectra of the p-Hj^O and H^^^O transitions (band 4b). Green 
lines show the model for SO2 emission from the hot core de- 
scribed in the text. 



in Figure 1, includes two relatively strong SO2 lines within the 
0-H2O spectrum, and a weaker line within the P-H2O spectrum. 
These lines are included, as an additional background, in the cal- 
culation of the H2O optical depth. Another SO2 line is present 
in the o-Hj^O spectrum at -15 kms ' (Fig. 2, upper panel) and 
has been similarly modelled. The full source model of S.-L. Qin 
(private comm.) suggests that no other strong lines are present in 
our spectra at velocities coiTesponding to the foreground clouds. 

Figure 2 (upper and middle panels) shows the o- and p-H^^O 
and Hj^^O spectra divided by the background emission, includ- 
ing dust continuum and the SO2 lines. The o-Hj^O is an equally- 
weighted average of the band la and lb data. The H^^O Unes are 
saturated over a wide range of velocities and thus unusable for 
a quantitative analysis. However, we have identified several ve- 
locity ranges with moderate saturation levels, marked with thick 
horizontal lines in Figure 2. These can be identified with the ex- 
panding molecular ring (< -50 km s~'), a transition between the 
4 kpc arm and the Orion arm (-12 to -7 kms"'), the Sagittarius 
arm (5 to 20 kms"'), and the Scutum arm (27 to 35 kms"'; pos- 
sibly blended with the Sagittarius B2 envelope; see e.g. Neufeld 
et al. 2000). Assuming that the foreground absorption com- 
pletely covers the continuum and all water molecules are in 
the ground state (a reasonable assumption for the diffuse fore- 
ground clouds given the very high spontaneous emission rate 
coefficients for the ground-state water lines, 3.458 x 10"^ and 
1 .842x 10"^ s"' for the ortho and para lines, respectively), we de- 
rive optical depths of the o- and p-water lines (t = - In I/Ig)- The 
resulting optical depth ratio is shown in Figure 2 (lower panel; 
left intensity scale). An ortho/para optical depth ratio of 1 corre- 
sponds to a column density ratio of 2. The resulting ortho/para 
column density ratio is given by the right hand scale in Figure 2 
(lower panel). 

The uncertainty in the line optical depth in a given velocity 
channel is given by 6t - exp{T)x6I/Io- The errorbars in Figure 2 
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Fig. 2. (Upper) Spectra of the ground state o- and p-H^^O lines 
towards Sagittarius B2(M) normalized by the continuum (green 
and red histograms, respectively). (Middle) Corresponding spec- 
tra of the Hj^O lines. (Lower) The ortho/para optical depth ratio 
(left scale) and column density ratio (right scale) as a function 
of LSR velocity. Black and blue points correspond to the H^^O 
and HPO measurements, respectively. 




Fig. 3. (Upper) Spectra of the ground state o- and p-Hj^O Unes 
towards W3 IC normalized by the continuum (green and red his- 
tograms, respectively). (Lower) The ortho/para optical depth ra- 
tio (left scale) and column density ratio (right scale) as a function 
of LSR velocity. 

Table 1. Column densities (cm-^) and ortho/para ratios towards 
Sagittarius B2(M). 



V(kms-') 


0/P 


MHoO) 


A'CHz) 


X(H20) 


-135 to -110 

-73 to -52 

-12 to -7 

5 to 20 

27 to 35 


2.35 ± 0.3 
2.35 ± 0.4 
2.8 ±0.5 
3.0 ±0.6 
2.3 ±0.3 


2.0 X 10'* 
3.7 X 10'* 
1.2x10'* 
6.7 x 10'5 
1.7x10'* 


5 X 10'" 

6 X 10'" 
4 X 10'" 
3 X 10" 
6 X 10'" 


4x 10-' 
6 X 10-' 
3 X 10-' 

2 X 10-* 

3 X 10-' 



(lower panel), are computed under the conservative assumption 
61 /lo - 0.05 (maximum uncertainty, dominated by the stand- 
ing waves, based on a comparison of the continuum normalized 
spectra of the 0-H2O line in bands la and lb; the signal-to-noise 
ratio for the P-H2O line is comparable, given the stronger con- 
tinuum, and the amplitude of the standing waves also scales with 
the continuum strength). Relative uncertainties of the o- and p- 
water optical depths have been added in quadrature. 

Table 1 gives weighted averages of the ortho/para col- 
umn density ratio in different velocity ranges towards 
Sagittarius B2(M), along with the coiTesponding uncertainties. 
Since the individual measurements are not independent but dom- 
inated by instrumental systematics, we use a more conservative 
estimate of the uncertainty of the mean for the two velocity inter- 
vals corresponding to the expanding molecular ring, computed 
from the peak-peak variation between the individual data points. 
The H2O (orthoH-para) column densities are also included, as- 
suming an H^^O/H^**0 ratio of 500 for the 5-20 kms-' com- 
ponent. We compute H2 column densities in the foreground gas 
using the method employed in Lis et al. (2001), based on '^CO 
absorption data, assuming a CO abundance of 1 x lO""* and a 
^^CO/'^CO ratio of 60 in the local gas in the Sagittarius arm 
(5-20 kms-' velocity range) and 30 in the remaining velocity 
intervals. The resulting column densities should be accurate to 
within a factor of 2. The H2O abundance in the various com- 
ponents is generally consistent with that derived by Neufeld et 



al. (2000; 4 - 7 x 10-'). The derived H2O abundance in the 5- 
20 kms-' component, based on the Hj^O measurements, is a 
factor of 3-6 higher than that derived for the other components 
based on the Hj^O data. However, the H^^O optical depth in this 
velocity range is low leading to large uncertainties. The water 
column density and abundance in this velocity range would be 
a factor of 2 lower if the gas responsible for the absorption is 
located in the Galactic center region (H^^O/Hj**© ratio of 250) 
rather than in the Sagittarius arms. 

The ortho/para ratio in the two velocity ranges between -12 
and 20 kms-' is consistent with the high-temperature limit of 
3, within the uncertainties. The ratio in the Scutum arm, at 27 
to 35 kms-' is lower than 3. However, this component may be 
blended with the Sagittarius B2 envelope. The low ortho/para ra- 
tios derived at the envelope velocities, from both H^^O and H^^O 
data, are likely caused by the excitation effects. Given the higher 
density of the gas, the assumption that all water molecules are 
in the ground state is no longer correct. The column density of 
0-H2O is more strongly affected, resulting in a lower apparent 
ortho/para ratio. In fact, we do see wings of the o-H^**0 line 
in emission at the envelope velocities (Fig. 2; middle panel). 
Deriving the ortho/para ratio in water in the Sagittarius B2 enve- 
lope requires detailed radiative transfer modelUng and is beyond 
the scope of the present Letter. 

We derive a low ortho/para ratio, 2.35 + 0.35, at velocities 
< -50 kms', corresponding to the expanding molecular ring. 
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In this case, the measurement uncertainties are small enough that 
the ortho/para ratio of 3 appears to be ruled out by our data. 

The o- and p- water spectra towards W31C, normalized by 
the continuum, are shown in Figure 3 (upper panel). Weak p- 
Hj^O absorption is seen at the cloud systemic velocity, but no o- 
or p-Hj**© absorption is detected in the foreground clouds. Once 
again, we have identified 3 velocity ranges where the o- and p- 
Hj^O lines are not completely saturated (thick horizontal lines in 
Fig. 3, upper panel). The resulting ortho/para ratio is relatively 
uniform, 2.8 + 0.2, consistent with the high-temperature ratio of 
3 within the measurement uncertainty (Fig. 3, lower panel). 



4. Discussion 

In the high-temperature limit one might expect an ortho/para ra- 
tio close to 3, when water is first formed. In the gas-phase, the 
excess energy of the exothermic reactions (e.g., recombination 
of HsO^) could lead to spin equilibration. For water molecules 
desorbed from ices on grain surfaces, one might also expect the 
initial ortho/para ratio to be 3 — if there is enough energy to des- 
orb a water molecule, either thermally or non-thermally, there 
is likely enough energy to populate numerous rotational states 
of the gas-phase species, both ortho and para, preserving a ratio 
close to 3 independent of the ortho/para ratio in the ice. However, 
the excess energy of formation is shared with the surface and the 
water ortho/para ratio may rapidly equilibrate at the grain tem- 
perature (e.g., Limbach et al. 2006). 

Once water molecules are in the gas phase, collisions with 
atomic and molecular ions (H^ and H^) can lead to proton ex- 
changes and, over time, produce a gas-phase ortho/para ratio be- 
low 3 and potentially commensurate w ith the gas tempera ture 
(see the discussion of H2 conversion bv lFlower et al.ll2006h . As 
the temperature approaches zero, all water molecules will be in 
the ground para state if collisions with ions can efficiently ex- 
change the ortho and para states. If there are no ions present, 
and therefore no barrierless collisions to exchange the ortho and 
para states of water, the ortho/para ratio of 3 will remain, inde- 
pendent of the temperature. Since ions are present in the ISM, 
we can expect the water ortho/para ratio to be generally lower 
than 3, given sufficient time. The time scale of the ortho/para 
conversion depends on the gas density. In dense clouds, there is 
a greater likelihood that the ortho/para ratio equilibrates at the 
gas kinetic temperature, as long as the density of ions is suffi- 
cient. In warmer, diffuse clouds, equilibrium would favor higher 
ortho/para ratios, but at lower densities equilibrium is less likely. 
Assuming a gas density of 10"^ cm"-*, a fractional abundance of 
protonated ions of 10"^, and a rate coefficient of 10"*^ cm^^s"', 
we estimate the time scale of the ortho-para equilibration to be 
of order 3 x 10^ years. 

We derive a low ortho/para ratio of 2.35 ± 0.35, coiTespond- 
ing to a spin temperature of ~27 K, towards Sagittarius B2(M) 
at velocities corresponding to the expanding molecular ring. The 
low ortho/para ratio may suggest water formation on dust grains, 
with water molecule spin populations equilibrated at the dust 
temperature. Alternatively, the high-temperature ratio of 3 may 
have relaxed in the gas-phase to a value in line with the ki- 
netic temperature of these clouds — a gas/dust temperature of or- 
der 30 K is quite reasonable for these relatively diffuse clouds 
within ~200 pc from the Galactic center (e.g., Tieftrunk et al. 
1994). The water ortho/para ratio has now been measured in at- 
mospheres of 9 solar system comets. The derived spin temper- 
atures range between 23 and 34 K, with an average of ~29 K 
(J. Crovisier, priv. comm.), close to our measurement in the ex- 
panding molecular ring. 



The higher ortho/para ratio, consistent with the statistical ra- 
tio of 3, that we derive towards Sagittarius B2(M) in the velocity 
range -12 to 20 kms"' may indicate higher gas or dust tem- 
peratures (e.g., Gardner 1988). Schilke et al. (2010) also derive 
a higher H20^ spin temperature towards this source at positive 
velocities, as compared to the line-of-sight clouds at negative ve- 
locities. In addition, strong chloronium absorption has been de- 
tected tow ards Sagittarius B 2(S) at velocities between -10 and 
20 kms~' (iLis et al.Ll2OI0l) . Since the H2C1^ abundance is en- 
hanced in warm, UV-irradiated regions, this indicates the pres- 
ence of a warm, diff'use gas component in this velocity range. 

In summary, in the diffuse regions studied here, both water 
formation mechanisms (gas and solid state) can contribute. With 
sufficient time (or higher densities with subsequently faster gas- 
phase reactions) chemistry in the gas can drive the spin tempera- 
ture towards the gas temperature. In addition, it has been postu- 
lated (Hollenbach et al. 2009) that photodesorption of water ices 
in low-density environments can release frozen molecules with a 
spin temperature that coincides with that of the dust grains. If the 
gas and dust temperatures are not equal this can lead to a ratio 
that reflects a mixture of these effects. Modeling the ortho/para 
ratio evolution is beyond the scope of this Letter, but the theme 
is that a ratio below 3 reflects cold environments (gas or dust 
temperatures below about 50 K). 

With improved calibration and better understanding of the 
instrumental effects, more accurate determination of the water 
ortho/para ratio in these and other sources will be possible in the 
future. However, the present work clearly demonstrates the out- 
standing spectroscopic capabilities of HIFI for providing robust 
constraints for the physical conditions and chemistry of the ISM. 
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